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We update the extraction of \Vus\ from hadronic r decay data in light of recent 
BaBar and Belle results on the branching fractions of a number of important strange 
\ decay modes. A range of sum rule analyses is performed, particular attention being 

paid to those based on "non-spectral weights", developed previously to bring the 
slow convergence of the relevant integrated D = 2 OPE series under improved con- 
trol. Results from the various sum rules are shown to be in good agreement with 
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avenues for future study are briefly discussed. 
PACS numbers: 12.15.Hh,13.35.Dx,11.55.Hx 

I. INTRODUCTION 

Once B factory data is fully analyzed, the precision of the hadronic r decay deter- 



! mination of 



Vus\ is expected to match or exceed that of K^^- and r[_ft'^2]/r[7r^2]-based 
methods [l|, 1], providing important tests of three-family unitarity and the Standard 
Model (SM) scenario for weak semi-leptonic AS* = 1 interactions. The analysis involves 
weighted integrals of the spectral functions of the vector (V) and/or axial vector (A) 
current- current correlators, examples of which are the normalized branching function ra- 
tios Rv/A;ij = r[r^ — > UT-hadionsv/ A;ij (7)]/r[T^ — > z/T-e~z/e(7)] associated with inclusive 
flavor ij = ud,us, V/A-current- induced hadronic r decays. 

The ratios Rv/A;ij are related to the J = 0, 1 spectral functions, Py/j^-ijis), of the 
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corresponding current- current correlators via 0], 
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where ijr = s/m^, wj.' {y) = (1 — + 2y), w}^' \y) = 2y{l — y)"^, Vij is the flavor 



ij CKM matrix element, and Sew is a short-distance electroweak correction. (0 + 1) 
denotes the sum of J = and 1 contributions. 

Using the general finite energy sum rule (FESR) relation, 
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valid for any correlator without kinematic singularities and any analytic weight, w{s), 

involving the 

correspondingly weighted linear combination of Ilyj^2,.j{s), the J 



the RHS in Eq. ([T]) can be replaced by an integral around \s\ = sq = 

0, 1 scalar parts of the 

relevant current- current correlator. For large enough sq, this integral can be evaluated 

<m2, either of 



using the OPE Similar FESRs can be constructed for any sq 



or sllyj^.-j{s) and general analytic w{s 



The corresponding spectral integrals will be 
denoted generically by i?J^(so) in what follows. We will also use the term "longitudinal" 
to refer to the J = contribution in any J = + 1/J = decomposition. 
Flavor-breaking differences of the ij = ud, us spectral data. 
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5]) provides a determination 



have OPE representations, 6Rqp^{so), beginning with a D 
ing \ Vud\ and any OPE parameters from other sources, Eq. 
of \Vus\ for any w{s), and any sq < sufficiently large that (i) convergence of the 
integrated series of known OPE terms is good and (ii) contributions from insufficiently 
well-known higher D OPE contributions are small. Explicitly, 
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For the weights used in the literature, at scales sq ~ 2 — 3 GeV^, SRqp^{so) is much 
smaller than Rudusi^o) (at the few-to-several- % level). The fractional uncertainty on \ Vus\ 
(~ A {6Rqp^{so)) /2R^^{sq)) produced by an uncertainty A ((^-RopbI'^o)) in ^RoPEi^o) 
is thus much smaller than that on 5R^p^{sq) itself, allowing modest precision in the 
determination of 5Rqp^{so) to be turned into a high precision determination of |V^s| [l|. 
Two important theoretical complications, however, must be dealt with. 

The first complication is the very bad convergence of the integrated longitudinal D = 2 
OPE series [HI and strong spectral-positivity-constraint violations produced by all trun- 
cation schemes employed in the literature for this badly converging series |6|. These 
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problems necessitate removing longitudinal contributions from the experimental spectral 
distributions so that one can work exclusively with sum rules based on the better-behaved 
J = + 1 OPE combination. Fortunately, for a combination of chiral and kinematic rea- 
sons, the spectral subtraction is strongly dominated by the well-known tt and K pole 
contributions. The remaining non-pole, or "continuum", contributions are doubly chi- 
rally suppressed, by factors of O ([m^ ± m^]^) (respectively, O ([m^ ± m„]^)) for the flavor 
ud (respectively us) V/A channels. The ij = ud non-pole longitudinal subtraction is thus 
numerically negligible. For ij = us, the continuum part of the longitudinal subtraction, 
though small, is not entirely negligible, but can be determined using input from sum 
rule and dispersive studies of the us scalar and pseudoscalar channels 0, [8|. The nor- 
malization of these results is strongly constrained by independent information on 
from the lattice and hence under good control theoretically. The resulting us continuum 
longitudinal subtraction is, as expected, numerically small. 

The second complication is the less than optimal convergence of the J = + 1D = 2 

OPE series. Defining An(s) = '^v+A-udi^) ~ '^v+A-usi^) ' exphcitly, with 

as{Q'^) and ms(Q^) the running coupling and strange quark mass in the MS scheme, 
and a = as{Q'^)/'K |9|, 

,2x1 OPE 3 ms{Q^^ ^ 



D=2 27r2 Q2 



1 + + 19.93a2 + 208.75a^ 



+ (2378 ± 200)a^ + ■■■], (5) 



where the coefficient is an estimate obtained using methods which previously provided 
a rather accurate estimate of the d^ coefficient in advance of the actual calculation of its 
value j^. Since independent high-scale determinations of as{Mz) correspond to a(m^) ~ 
0.10, convergence at the spacelike point on \s\ = sq < rn^ is marginal at best. While the 
decrease of |a<j(Q^) | on moving away from the spacelike point along the contour does allow 
the convergence of the integrated series to be improved through judicious choices of the 
sum rule weight w(s), for weights not chosen specifically with this purpose in mind very 
slow D = 2 convergence is to be expected. Such slow convergence is, for example, seen 
explicitly in the case of the so-called {k,m) spectral weights t9|. Because of the growth 
of as with decreasing scale, premature truncation of such a slowly converging series will 
produce \Vus\ values with an unphysical so-dependence. In view of the intrinsically slow 
convergence of the underlying D = 2 correlator series, a study of the so-dependence of 
the output \Vus\ (and demonstration of the existence of an so-stability window) is thus 
crucial to establishing the reliability of any hadronic-r-decay-based FESR extraction of 
\Vus\- Analyses in which such a study has not been performed thus typically contain 
additional theoretical systematic uncertainties, which may be significantly larger than 
the theoretical errors explicitly studied. 
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TABLE I: ^us^r branching fractions for strange hadronic states X~g. Column 2 contains 

the PDG2006 world averages, column 3 the 2007 Babar and/or Belle results, and column 4 the 
resulting new world averages, with PDG-style S factors included where needed. 
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II. THE UPDATED DETERMINATION OF \Vus\ 

Our results for \Vus\ are based on the ALEPH ud [lol] and us spectral data, for 
which information on the relevant covariance matrices is publicly available. Since no 
V/A separation has been performed for the us data, we work with the V+A combination 
(the ud spectral integral errors are also reduced by working with the V±A combination). 
The older us data has been modified to incorporate the recent (2007) BaBar 12 



and Belle [13| updates of a number of the strange branching fractions, Bx^^- The latter 
measurements, together with the previous and resulting world averages, are shown in 
Table [B Since the re-measurement of the strange spectral distribution is not yet complete. 



we follow Ref. [ij] and obtain a partially updated inclusive sum by rescaling the 1999 



ALEPH distribution [11[, mode- by-mode, by the ratio of new to old B values. 

With the above ud and us V±A spectral distributions as input, \Vus\ is determined, as 
a function of the upper integration limit, sq, in Eq. ([2]), using J = ± 1 FESRs based on 
(i) the three improved-convergence, non-spectral weights, wio, wio and W20, introduced in 
Ref. 15], and (ii) the kinematic weight (expected to display problematic integrated 



OPE convergence behavior and hence so-instability for \ Vus\ [l, l^)- Standard OPE input, 
the correlator form of the D = 2 contributions, and input for the K pole contribution 
to the us spectral integrals are employed. More details of the analysis will be reported 



elsewhere [17[. The resulting sq dependence of \Vus\ is shown in Fig. [H 
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The figure shows reasonable So-stabihty for wiq, Wiq and W20 in the region above 
So ~ 2.5 GeV^, but significant so-instabihty for w!^'^\ The Sq = |\4s| values, 

0.2156 ± 0.0028exp ± 0.0022*,, for W20 

0.2154 ± 0.0032e^p ± 0.0015*;, for ww 

0.2149 ± 0.0033exp ± 0.0010*,, for wio 

0.2144 ± O.OOSOexp ± 0.0017*,, for wj*'"^ , (6) 



are in good agreement, but lie ~ 3a below the expectations of 3-family unitarity |1( 
The theoretical errors have been estimated using the scheme discussed in Refs. [2I , which 
is, by design, more conservative than that employed in previous analyses. In spite of this, 
the So-instability of the results from the w!^'^^ analysis suggests that the theoretical error 
in this case may still be insufficiently conservative. The ~ 0.0050 downward shifts in 
\Vus\, compared to the 2006 versions of these analyses [H, [2|, result dominantly from the 
-0.075% (-0.050%) shifts in Bktt {Bx-n+n-), which represent, respectively, 2.6% (1.8%) 
decreases in the total strange branching fraction and hence, from Eq. (jlj), ~ 1.3% (0.9%) 
downward shifts in \Vus\- 

In Figure [2] we compare the results above for \Vus\ to those of a number of other 
determinations, from various sources. The result from Ki^ decay is based on the most 
recent KLOE K decay data 3], together with the RBC/UKQCD lattice determination 



of /+(0) [19^], that for Marciano's method (which employs the ratio of Kfj_2 to 77^2 decay 
widths in combination with lattice input for fxl fir |20|) on the FlaviaNet Kaon Work- 
ing Group assessment of the experimental ratio |2l| and Juttner's Lattice 2007 plenary 
assessment of the status of lattice determinations of fKlfv |22|, while that from hyperon 
decays, included for completeness, is Jamin's assessment as of Electroweak Moriond 2007 
(which is 0.0010 higher, and with larger errors, than the result quoted in Ref. [23] )• The 
first of the r decay results represents a range of assessments based on the FESR 
and pre-2007 us data. The next two FESR assessments are from Ref. [23] and 

incorporate most, but not all, of the new 2007 us data. The sizeable difference produced 
by switching from K^2 input for the K pole contribution (the "pred. r — > Kz/" result) 
to the central measured value (the "meas. r Ku" result) is a reflection of the high 
degree of cancellation between the ud and us spectral integrals [2^. Note that the r 
decay results just mentioned, as well as those of Refs. [26|, employ the 4-loop-truncated 
Adler function evaluation of the D = 2 OPE integral contribution. The remaining r 
decay results shown in the flgure are the ones obtained by us above, which, in contrast, 
employ the 4-loop-truncated correlator evaluation of the D = 2 contribution. The non- 
trivial difference between the results obtained using the two different methods (which are 
formally equivalent, up to contributions of O(a^) and higher) provides further evidence 
for the slow convergence of the integrated D = 2 series for the wl^'^^ weight, seen already 
in Fig. [TJ 

While the results above are compatible with small beyond-the-SM contributions, such 
a conclusion is obviously premature until BaBar and Belle have cross-checked each other's 
recent results and also flUed in the missing entries in column 3 of Table [T] (especially for 
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FIG. 1: \Vus\ versus sq for, from top to bottom, ii;20i ^^lO) ^^lo and wj.' . 

the sizable ^°7r°7r~ and previously estimated, but unmeasured, KSn and KAtt modes). 
Results for higher multiplicity modes, with branching fractions down to the few xlO~^ 
level, are also needed for a determination at the desired level of accuracy. A finalized 
version of the results will require completion of the re-measurement of the full us spectral 
distribution, work on which is in progress. Note that improved precision on the measure- 
ment of the ratio of branching fractions of the single-particle vr and K modes is also of 
interest for the determination of \Vus\ since it would allow an independent detemination 
in the r sector analogous to that based on r[i^^^2]/r[7r^2] [s^]- Experimental input on the 
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FIG. 2: \Vus\ from various sources 



chirally-suppressed continuum us J = spectral contributions, needed for performing 
the longitudinal subtraction, appears potentially experimentally feasible, and would also 
be most welcome. 

Finally, we note that a re-measurement of the ud distribution is also of relevance, par- 
ticularly in light of the known discrepancy between the tttt and 47r distributions measured 
in T decay and those implied by electroproduction cross-sections (via CVC and known 
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isospin-breaking corrections). Were the electroproduction expectations to be correct, 
the values of \Vus\ obtained from the above analyses would be increased by ~ 0.0018. 
It is worth noting, in this regard, that two recent analyses of non-strange hadronic r 
decay data [27|, i28|], one designed specifically to incorporate an improved treatment of 



D > A OPE contributions which were problematic for earlier determinations, yield val- 
ues of as{Mz) in excellent agreement with the results, 0.1190(26) and 0.1191(27), ob- 
tained frona updates to the global fit to electroweak observables at the Z scale (quoted 
respectively). Two recent updates jsil, 3^ of the earlier lattice determi- 
employing slightly different implementations of the original approach, also 
obtain values (0.1183(7) and 0.1192(11), respectively) which are higher than the original 
result and now in excellent agreement with the results of the r extraction and global 
electroweak fit. Since electroproduction cross-sections imply a significantly lower value. 



in Refs. 
nation 13 




as{Mz) ~ 0.1140 — 0.1150 [34], we consider it unlikely that the ud r data is significantly 



in error, but nonetheless a re-measurement from the B factory experiments should be 
pursued, given the potential impact on both SM expectations for {g — 2)^ and on the 
value of \Vus\ obtained in the hadronic r decay analyses. 
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